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STMIARY 

We  oonslder  Earth  satellites  in  the  region  where  the  perturbing  effeots 
due  to  Earth's  ob lateness  and  luni-solar  gravitational  forces  are  oonparable. 

A  general  solution  is  obtained  for  simultaneous  preoession  about  any  mnber  of 
fixed  axes;  this  is  an  extension  of  Laplaoe's  treatment  far  the  motion  of 
Iapetus  about  Saturn.  Results  arc  given  for  general  orbits  on  the  assumption 
that  the  lunar  orbit  lies  in  the  eoliptio.  Synchronous  orbits  are  considered 
in  greater  detail. 
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1  INTRODUCTION 

The  main  perturbations  for  a  relatively  distant  Earth  satellite  are  due 
to  the  Barth's  oblate  ness  and  the  gravitational  attractions  of  the  Sun  and  Boon. 
We  consider  satellites  at  distanoos  of  a  few  Earth  radii,  where  the  effeot  of 
the  Earth's  oblate  ness  is  of  the  same  order  as  the  effeota  of  the  distant 
bodies.  This  range  oovers  a  major  region  of  interest  for  oaomunioation 
satellites.  Two  aspeots  of  this  problem  have  alroady  been  considered  by  the 
authors  separately1  >2,  and  the  present  work  gives  a  more  ocmpleto  solution  of 
the  general  problem. 

The  olassioal  problem  of  this  type  oonoerns  the  motion  of  Iapetus^  which, 
apart  from  the  very  distant  and  retrograde  satollite  Phoebe,  is  the  outermost 
satellite  of  Saturn.  The  perturbing  force  due  to  tho  ob lateness  of  Saturn 
(suitably  modified  to  inolude  the  attraction  of  the  rings  and  the  inner 
satellites)  is  oomparable  with  the  perturbing  foroe  due  to  tho  gravitational 
attraction  of  the  Sun,  By  itself  the  former  would  cause  tho  satellite's  orbital 
plane  to  process  around  the  axis  of  rotation  of  Saturn,  while  the  latter  by 
itself  would  oause  the  orbital  plane  to  preoess  around  the  pole  of  the  orbit  of 
Saturn.  Tho  motion  of  Iapetus,  which  is  a  oonbination  of  simultaneous  prooea- 
sions  about  two  different  axes,  was  first  considered  by  Laplaoe^  and  later 
discussed  in  considerable  detail  by  Tisseran&5, 


We  take  the  disturbing  function  to  inolude  the  second  sonal  harmonic  of 
the  Earth's  gravitational  potential  and  the  lowest  order  term,  i. e.  the  seoond 
harmonic,  in  the  expansion  of  each  of  the  solar  and  lunar  disturbing  functions. 
We  consider  only  initially  oiroular  satellite  orbits  and,  since  the  seoond 
harmonic  has  no  secular  effeot  on  the  eooontrioity  of  a  oiroular  orbit,  the 
orbits  will  remain  oiroular  throughout.  In  faot  the  higher  harmonies  in  the 
lunar  disturbing  function  become  more  important  as  the  orbital  radius  increases, 
and  their  neglect  sets  an  upper  limit  to  the  radius  of  the  satellite  orbit  for 
which  the  theory  is  applicable,  especially  since  the  parallaotio  term  (third 
harmonio)  will  ohange  the  eooontrioity  of  a  oiroular  orbit.  A  small  initial 
eooontrioity.  such  as  may  be  built  tg>  by  the  parallaotio  term,  may  increase 
rapidly  if  the  orbit  is  sufficiently  large6. 

The  disturbing  functions  aro  averaged  analytically  over  the  mean  anomalies 
of  both  the  satollite  and  the  disturbing  bodies.  Par  this  procedure  to  give  a 
good  representation  of  the  satellite's  motion,  the  period  of  the  predioted  long- 
period  motion  of  the  orbital  plane  should  be  reasonably  long  compared  with  the 
periods  of  the  disturbing  bodies. 


As  a  result,  of  tho  limitations  discussed  in  the  last  two  paragraphs  our 
treatment  is  most  useful  for  satellites  at  distances  between  about  3  and  10 
Earth  radii. 


DgTQBBire  WTNOTIQN  AND  EqOATIONS  CP  MttKW 


The  perturbing  effeot  on  the  motion  of  a  satellite  about  the  Barth  due  to 
the  gravitational  attraction  of  a  third  body  of  mass  m^  ( j  ■  1  for  Sun,  j  ■  2 

for  Moon)  is  given  by  the  disturbing  function 


(D 
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where  r,  r ^  are  the  poeitions  of  the  satellite  and  the  disturbing  body  relative 
to  the  Earth,  ani  (i .  ■  Gm..  We  assvsne  r  <<  r.  and  retain  only  the  lowest -order 

0  0  0  j 

term  in  the  expansion  of  this  disturbing  function  in  powers  of  r/Tjf  via* 


P2(ooe 


.22 

•irrj 


(2) 


where  0^  is  the  angle  between  r  and  r^,  Including  the  seoond  harmonic  in  the 
Earth's  gravitational  potential,  the  oonplete  disturbing  function  may  be  written 


Here  ^  •  Gm,  where  m  is  the  mass  of  the  Earth,  Jg  is  the  ooeffioient  of  the 
seoond  sonal  hannonio,  Rg  is  the  mean  equatorial  radius  of  the  Earth,  and  6q 
is  the  angle  between  r  and  the  polar  axis  of  the  Earth. 

Par  the  long-period  motion  it  is  permissible  to  average  the  disturbing 
function  over  the  mean  anomalies  of  the  satellite  and  of  the  disturbing  bodies. 
Equation  (3)  Is  first  written  in  the  form 


0 


2  2 


]> 


(4) 


where  R  is  a  unit  veotor  along  the  Earth's  axis.  To  average  over  the  mean 

O  mm 

anomaly  l  of  the  satellite  the  follcsring  integrals  will  be  used  (of.  Ifusen'). 


‘  •  (5*> 

2x 

&f  ^5-£d4  "  &-&B)  ,  (5b) 
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&  i  r2«  ■ 


2% 

h  J  zzu  a 


where  a,  e  are  the  semi-major  axis  and  eooentrioity  of  the  satellite's  orbit, 
and  P,  £,  R  are  the  natural  unit  vectors  associated  with  the  satellite  orbit, 
with  R  along  the  positive  normal  to  the  orbital  plane,  P  in  the  orbital  plane 
in  the  direction  of  perigee,  and  £  =  R  «  P.  is  the  iSemtensor.  Equations 

(5a,b)  may  easily  be  verified  taking  the  true  anomaly  as  the  variable  of  inte¬ 
gration,  and  equations  (6a,b)  by  taking  the  eocentrio  anomaly.  It  is  immaterial 
in  what  order  the  double  averaging  is  performed,  and  it  proves  rather  simpler  to 
average  first  over  the  mean  anomalies  of  the  disturbing  bodies  using  the 
analogues  of  (5s)  and  (5b)  so  that  the  singly-averaged  disturbing  function  may 
be  written  as 


a2(l+i  e2^  ,  (6a) 

*  a2  f(l  +  l*e2)  PP  +  (1  -  e2)  ,  (6b) 


U 


“j(1  - 


(7) 


where  a^,  e ^  are  the  semi-major  axis  and  eooentrioity  of  the  orbit  of  the  j**1 
disturbing  body,  and  is  the  unit  vector  normal  to  the  orbital  plane,  i.  e. 

R^  lies  along  the  pole  of  the  eoliptio  and  Rg  is  normal  to  the  lunar  orbital 
plane. 

After  taking  the  average  over  the  mean  anomaly  of  the  satellite  using  (5a) 
to  (6b),  the  doubly-averaged  disturbing  function  takes  the  form 


U  »  n  a  U*  , 


(8) 


where 


U* 


«„<’  -  •2>"3A  (ks  .  ^>‘  -  f] 

*  T  “i  [«'  -  *2)(s  •  »/  ♦  «2  [i  - 1  (£  .  a/]]  ,  (9) 


where  n  is  the  mean  motion  of  the  satellite,  and 

"  5  “ 
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(10) 


The  aeni-toaijor  axle  a  and  the  mean  motion  n  are  oonstant  in  the  long-period 
motion,  end  they  have  been  freely  incorporated  in  (10),  In  writing  (9)  we  have 
used  the  tensor  identity 


PP+23+RR  >  ^ , 


and  dropped  an  irrelevant  oonstant  tern  from  the  luni-solar  part  of  the 
disturbing  function. 

If  we  ohoose  a  particular  reference  plane,  we  oan  express  the  disturbing 
funotion  in  terms  of  the  corresponding  elements  e,  i,  0,  u  and  give  the  Lagrange 
equations  for  their  variation.  However  this  prooedure  is  necessarily  very 
complicated  if  the  two  disturbing  orbit  planes  are  taken  as  distinct,  and  it  is 
muoh  simpler  to  use  vsotorial  elements  whioh  are  not  tied  to  any  particular 
frame.  The  veotorial  elements  are  taken  as  e  and  h  defined  by 


e_  ■  e  P  ,  and  h 


(1  -  e2)*R 


(11) 


so  that 


0,  and  e  ♦  & 


(12) 


In  terms  of  these  elements  the  Lagrange  equations  take  the  ooopaot  and 
symnetrioal  form® 


h 


(13) 


where  U*  a  U/ne^s  From  (9)  the  disturbing  function  U*  may  be  written  in  tenet  of 
these  elements  as 

r  o  "5/2  9  a  o  -3/2 

U*  -  W0  j*0  -  e2)  (£ .  r/  -  i  (1  -  e2)  J 

♦  ^  [Wfc .  Sj)2  ♦  «2  -  I  (£  .  S,j)2]  .  (14) 

>1,2 
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Prom  (13)  and  (14)  the  equation  for  h  is 

.  2  "5/2 

h  =  -  «0(1  -  ez)  (h  .  Rq)^  «  h) 


• )  ‘  “  j(h  •  Lj)(R.j  -  h)  +  5  V  «j(e  ,  RpCRj  -  e) 

£i'f2  M2 


(15) 


It  ia  obvious  from  the  form  of  U*  that  dU*/d£  vanichea  for  e  a  0,  so  that  £ 
must  vanish  for  a  oiroular  orbit  (with  our  assumptions  regarding  the  disturbing 
function). 

We  near  consider  only  circular  satellite  orbits,  setting  e  a  0  in  (15)  and 
replacing  h  by  R,  so  that  we  can  write 

2 

R  =  -  Y  «j(R  .  Lj)(Rj  -  R)  .  (16) 


If  only 
normal  to  the 

constant  rate 

precession  is 


one  of  the  terms  on  the  right  hand  side  of  (l6)  were  present f  the 
orbit  R  would  simply  regress  around  the  appropriate  R^ata 

u  j  00s  Yj,  where  y j  ia  the  angle  between  R  and  R^.  The  rate  of 

zero  when  Y  j  ■  90°. 


For  a  distant  Earth  satellite  the  motion  of  the  orbital  plane  is  a  com¬ 
bination  of  rates  of  procession  about  three  differont  axes.  Moreover  the  lunar 
orbital  plane  ie  inclined  to  the  eoliptio  at  5°9'  and  is  itself  regressing 
about  the  pole  of  the  eoliptio  with  a  period  of  18.6  years.  It  is  usually 
adequate  to  approximate  by  assuming  that  the  lunar  orbit  lies  in  the  eoliptio, 
or  alternatively  by  taking  a  mean  position  for  the  pole  of  the  lunar  orbit  if 
the  time  interval  of  interest  is  short. 


3  GENERAL  SOLUTION  FOR  AMY  NUMBER  OF  FIXED  AXES 


We  shall  show  that  equation  (16)  can  be  solved  exactly  for  simultaneous 
preoession  about  any  number  of  fixed  axes.  However  we  first  oonsider  the  simple 
approximation  of  as  aiming  that  the  angles  y,  between  R  and  the  various  axes  R . 

V  V 

are  all  small.  This  requires  that  the  axes  Rj  themselves  must  be  dose  together, 
which  is  approximately  true  in  practioe  since  and  R^  are  inclined  at  23°, 

R^  and  Rg  at  5°,  and  R^  and  Rg  at  between  18°  and  28°,  If  the  various  ooslne 
factors  (R  .  R^)  »  00s  y j  are  simply  replaced  by  unity,  equation  (16)  becomes 


it* 


55a, 


(17) 
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where 


u  R 


udh  ’ 


(18) 


and  R  ia  a  unit  vector.  Then  R  regresses_at  a  oonstant  inclination  around  the 
mean~pole  R  with  the  preoession  velocity  «. 

To  obtain  an  exaot  solution  of  (l6),  we  note  that  the  equation  simplifies 
considerably  if  it  is  referred  to  the  principal  axes  of  the  synmetrio  tensor 
£  « j  Rj  Rj  .  Alternatively  one  can  ask  whether  there  are  any  positions  whore  R 

will  remain  at  rest.  In  fact  ft  will  vanish  if  tho  scalar  product  of  R  with  the 
tensor  either  vanishos  or  is  a*multiple  of  R,  In  either  oaso  R  is  an~oigen- 
veotor  of  the  tensor,  i, o.  ~  ~ 


(l  Mj  Rj  Rj  *  X  j)  ,  R  a  0  , 


(19) 


where  1_  is  the  idem  tons  or,  and  X  is  tho  eigenvalue. 

To  obtain  an  idea  of  the  magnitudes  involved,  the  eigenvalues  have  beon 
evaluated  using  the  approximation  that  the  lunar  orbit  lies  in  the  eoliptio, 
i.e.  that  Rg  coincides  with  R^ .  The  Rj  are  first  referred  to  a  right-handed 

system  of  equatorial  axes  with  one  axis  towards  the  vernal  equinox.  The 
oenponents  of  are  (0,0, i),  while  those  of  R^  are  (0,  -sin  c,  oos  e),  e  being 

the  obliquity  of  the  eoliptio.  Then,  if  we  set  +  «2  »  u*,  the  eigenvalues 
are  the  roots  of  tho  equation 


—X 


«*  sin  e  -  X 


sin  e  oos  e 


0  -«*  sin  e  oos  e  w  +  «*  oos  e  -  X  ! 

o 


o  . 


We  label  the  eigenvalues  X^ ,  Xg  and  X  ^  and  arrange  them  in  the  order 

Xj  >  Xg  >  X^,  The  suffixes  1,  2  and  3  used  here  for  the  eigenvalues,  and  later 

for  the  principal  axes,  have  no  connexion  with  the  suffixes  1  and  2  used  for  the 
Sun  and  Moon,  We  find 

N  »  o. 

XjJ  ■  i(wo  +  «*)  *  i^“0  +  *»*)2  -  W*  Sin2cj^  , 

The  values  of  «0,  «1  and  «2  and  the  eigenvalues  Xg  and  Xj  are  plotted  as  functions 
of  a/^  in  Rig.  1. 
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The  directions  of  the  corresponding  principal  axes  are  readily  discovered 
from  (16),  That  oorreeponding  to  the  eigenvalue  aero  is  olearly  perpendicular 
to  R  aid  1^,  i»e.  it  is  direoted  towards  (or  away  from)  the  vernal  equinox. 

The  remaining  two  principal  axes  lie  in  the  plane  of  Rg  and  R^.  If  «  is  the 
angle  between  R  aid  the  principal  axis  corresponding  to  (or  A,,),  we  have 


tan  2«  ■  u*  sin  2s/(uq  +  u*  oos  2«)  . 


The  orientation  of  the  principal  axes  0123  is  shown  incidentally  in  Fig.3»  and 
the  value  of  the  angle  a  is  given  as  a  function  of  a/Rj,  in  Fig.  2. 

Assuming  that,  in  the  general  oase,  we  have  found  the  principal  axes  01, 
02,  03  aid  the  eigenvalues  A^,  Ag,  Aj,  equation  (l6)  beoomes  in  soalar  form 

X1  *  (X3  "  V  *2  x3  I 

| 

*2  »  (X1  -  X3)  x3  x,,  l  (20 i 

*3  *  (X2  -  xi)  xi  *2  )  » 

where  x^ ,  Xg,  Xj  are  the  components  of  R  along  the  principal  axes.  Multiplying 
each  equation  by  the  oorreeponding  x^  and  adding,  we  have 

*1  *1  *  *2  *2  +  *3  *3  *  °  * 

which  merely  leads  to  the  known  condition  that  R  is  a  unit  vector: 


+ 


■f  X, 


(21) 


On  the  other  hand  multiplying  each  equation  by  X^  x^  and  adding  gives 
X1  xl  X1  ♦  h  *2  *2  +  X3  x3  x3  “  0  * 


which  yields  the  non- trivial  integral 

X1  X1  +  X2  *2  +  X3  *3  "  Xo 


(22) 


where  XQ  is  a  constant  lying  between  A,,  and  Xj,  from  (21 ). 
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The  integral  (22)  oan  he  found,  without  diagonalieing  the  tensor,  in  two 
other  ways.  First,  multiplying  (l6)  soalarly  by  ^  uk(R  ,  R^)  R^  gives 

k 


Y  *k®  •  4X4  •  V  ■  Y  “s  •  V®  •  4>  £  *  Sj>  •  4] 

k  J,k 


The  right-hand  side  of  this  equation  must  vanish  since  it  is  anti-syianetrlo  in 
j  and  k  so  that  we  obtain 


2 

ooe  Yj 


(23) 


whioh  is  the  form  in  which  the  integral  is  given  for  simultaneous  preoession 
about  two  different  axes  in  the  classical  problem  of  the  motion  of  Iapetus3, 


Seoondly,  we  note  that  the  doubly-averaged  Hamiltonian  of  the  system  is 
H  ■  -(|i/2a)  -  U,  If  the  axes  R^  are  fixed,  then  K  is  independent  of  explioit 

tine  and  must  be  a  oonstant  of  the  motion.  Sinoe  a  is  also  a  oonstant  of  the 
long-period  motion,  U  itself  must  be  constant.  Using  the  knowledge  that  an 
initially  oiroular  orbit  remains  circular,  we  recover  the  integral  (22)  on 
setting  e  •  0  in  (9). 


4  NATURE  OP  THE  MOTION 

The  equations  (20)  are  similar  to  the  Euler  equations  for  the  motion  of  a 
rigid  body  under  no  foroes,  and  corresponding  results  oan  be  derived  for  the 
nature  of  the  motion  and  the  period.  Ue  reoall  that  the  eigenvalues  are 
arranged  in  the  order  Xj  >  X2  >  X1  and  that  X0  lies  between  X1  and  X,,  In 

praotioe  X^  is  muoh  greater  than  the  other  two  eigenvalues  (see  Pig.1), 

Equation  (22)  gives  a  one-parameter  family  of  curves  on  the  unit  sphere  which 
are  the  possible  trajectories  of  the  pole  Rj  the  initial  position  of  the  pole 
determines  the  value  of  XQ,  The  motion  is  periodic,  and  tho  trajectories 

divide  into  two  types  depending  on  whether  \Q  is  greater  or  less  than 

Case  (l)i  If  Xj  >  XQ  >  Xg,  (21 )  and  (22)  cannot  be  satisfied  simultaneously 
if  x^  m  0,  so  that  the  pole  R  oan  never  attain  or  oross  the  plana  x^  ■  0, 

Case  (ii)t  Similarly  if  Xg  >  XQ  >  X^,  the  pole  oan  never  attain  or  oross  the 
plane  x^  ■  0. 
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Eliminating  Xg  botwcen  (21 )  and  (22) ,  we  find 

(Xj  -  Xg)3^  "  (^2  "  ^1  *  (*0  "  ^2^  9  (^) 

so  that,  on  setting  ^  the  bounding  curve  between  the  two  types  of 
trajectory  is  the  intersection  of  the  pair  of  planes 

-  t  J(a2  -  x1)/(x3  -  x2)J*  (25) 

with  the  unit  sphere.  The  half -angle  between  the  planes  ia  plotted  (under  the 
legend  x/2  -  $2)  in  fig, 2  as  a  function  of  a/Rg. 

An  idea  of  the  shapes  of  the  trajectories  is  given  by  their  projections 
on  the  coordinate  planes;  eliminating  *j  and  x^  in  turn  from  (21 )  and  (22)  we 

find 

(X^  *  X^  )x^  +  (Xj  ~  X2)x2  ■  (X^  -  Xj)  (26a) 

(X2  "  X^ )Xj,  ♦  (Xj  -  X^ )Xj  a  (x^  ■  X^)  ,  (2ft) 

for  Xj  >  XQ  >  Xg,  the  trajectory  encircles  the  axis  03,  and  its  projection  on 

012  is  an  ellipse  of  small  eccentricity  (assuming  X^  »  Xg,  X^)«  Likewise  for 

X2  >  X©  >  X^,  the  tra jootory  onoirolos  the  axis  01 ,  and  its  pro  jeotion  on  023 

is  an  ellipse  of  large  eooentrioity,  Trajootorios  of  the  pole  R  are  shown  in 
Pig.  3. 

Laplace  considered  the  oaso  where  the  polo  of  the  satellite  orbit  is  very 
near  the  prinoipal  axis  03,  when  the  trajectory  on  the  unit  sphere  is  very  nearly 
a  oirole  around  the  pole  03.  This  is  equivalent  to  Laplace's  result,  which  is 
that  the  satellite  orbit  regresses  at  nearly  oonstant  rate  and  inclination  on  the 
proper  or  Laplaoian  plane  which  is  here  simply  the  coordinate  plane  012, 

5  PERIOD  OP  qqr.TT.TJLTTnN 

Using  (26a)  and  (2ft)  to  eliminate  x,  and  x.  from  the  seoond  of  (20).  we 
find  1  5 

*2  3  1  ^3  “  \>)  ~  (^3  "  [Xe  **  X^)  “  (X2  “  X^)*^J  « 
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The  quartor-poriod  of  the  osoillation  is  given  by  inverting  and  integrating  from 
Xg  *  0  to  *2  ■  th0  uPPer  linit  as  a  function  of  Xg  depends  on  whether 

X  >  X,  or  not,  For  instanoe  if  >  X«,  the  uppor  limit  of  integration  is 

determined  by  tho  first  square  root,  and  the  period  of  a  complete  osoillation  is 

v^ 

T  .  k  T(X3  -  XQ)(X0  -  X,)]4  j  (1  -  x^/v2)"*  (1  “  **2  » 

o 

where  v2  >  (X^  -  XQ)/(X3  -  X2)  , 

aid  v2  .  (X0  -  X1)/(X2  -  X^  . 

The  transformation  «  ^  sin  h  reduces  tho  above  to  the  standard  form 

*/2 

T  a  V(X3  -  X2)"*  (X0  -X^f  (1  -  k2  sin2t)  a*  , 

o 

where 

k2  .  v2/v2  a  -  Xq)(x2  -  X^j/^  -  X2)(Xo  -  X^J  , 

so  that 

T  .  V(X3  -  X2)"^  (Xo  -  X^  K(k),  Xo  >  X2  ,  (27a) 

where  K(k)  is  the  oornplete  elliptic  integral  of  the  first  kind.  As  \Q 

approaches  its  upper  limit  the  trajectories  shrink  towards  the  stationary 

point  on  the  axis  03;  k  tends  to  sero,  and  the  period  tends  to  a  limiting  value 
which  we  will  call  T3>  given  by 

Tj  a  2x(X3  «•  X^  )"^  (Xj  ~  Xj,)"^  ,  (28a) 
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slnoe  K(0)  *  r/2.  On  the  other  hand  as  Xq  approaches  Xg,  k  tends  to  1,  and  the 
period  tends  to  infinity. 

The  period  can  he  evaluated  similarly  for  Xq  <  X2  giving 

T  »  4(X3  -  Xo)"*  (X2  -  X,)*  K(k'),  \  <  X2  ,  (27b) 

where 

k'2  a  |^(X3  -  X2)(Xq  -  N “  Xo^X2  "  * 

As  X0  approaches  its  lower  limit  X^ ,  the  trajectories  shrink  towards  the  stable 

point  on  the  axis  01;  k'  tends  to  zero  and  the  period  tends  to  the  limiting 
value  T^  where 

T1  «  2x(X3  -  X,,)^  (X2  -  Xt)“*  .  (28b) 


Pig,  4  gives  the  limiting  periods  and  T^  as  functions  of  a /teg,  again  on  the 
assumption  that  the  lunar  orbit  lies  in  the  eoliptio  (X^  a  0). 
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>U3  CRB  ITS 


For  orbits  with  a  period  of  one  sidereal  day  (at  a  height  a/kg  a  6. 6108) 
we  have  the  following  values 


wo  a  4. 900°/yr. ,  w1  a  0. 738°/yr, ,  «2  •  1.6l1°/yr. 

t/e  consider  first  the  approximation  that  the  lunar  orbit  lies  in  the 
eoliptio,  so  that  the  satellite  orbit  is  subjected  to  simultaneous  precession 
foroes  about  only  two  axes,  the  Earth's  axis  and  the  pole  of  the  eoliptio.  . 

If  we  take  the  simplest  approximation  of  equation  (17),  then  the  mean  pole  R 
lies  in  the  plane  of  R  and  R.  at  an  angle  of  7°33'  from  the  Earth's  axis,  and 

•  -O  ^  "  • 

u  a  7,11o /yr.  The  shortest  period,  when  R  is  very  olose  to  R,  is  50.6  yr.  If 
the  satellite  orbit  is  initially  equatorial,  the  rate  of  preoession  is  reduced 
by  the  oosine  of  the  inclination  to  fi,  and  the  period  is  51,0  yr. 

In  the  general  treatment  of  Section  3,  but  still  assuming  the  lunar  orbit 
lies  in  the  eoliptio,  the  eigenvalues  are 

X1  a  0,  X2  a  0,26l°/yr. ,  X3  ■  6, 988°/yr. , 
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aid  the  principal  axis  03  liee  in  the  plane  of  ^  and  at  an  angle  a  a  7°23' 

to  the  Barth's  axis.  Prom  (28a)  the  limiting  value  of  the  period  when  R  is  very 
olose  to  03  is  52.5  yr.  On  the  other  hand  if  the  satellite  orbit  is  initially 
equatorial, 


\Q  a  6.877°/yr.,  k2  a  6.26  x  10’3  , 

and  the  period  is  52, 9  yr.  Thus  the  simplest  approximation  of  a  mean  pole  is 
ren*rk&bly  accurate  for  this  oaae. 

The  period  of  the  osoillation  oan  be  oaloulated  from  (27s)  and  (27b)  for 
all  possible  trajectories  at  synchronous  height.  Since  all  trajectories 
intersect  the  plans  013,  a  convenient  parameter  here  is  the  angle  0  between  03 
and  the  intersection  of  the  trajectory  with  013,  Setting 

X1  a  Sin  Xg  a  0,  Xj  a  OOS  <f>  , 

in  (22)  we  have 

2  2 

X  a  X.  sin  6  ■ f  X,  oos  6  . 

o  1  r  3 


The  value  of  ^  corresponding  to  Xq  a  X^  is  denoted  by  ^  so  that  the  half-angle 

between  the  bounding  planes  between  the  two  regions  of  motion  in  Pig.  3  is  */2  -  ^ 

(which  has  already  been  plotted  in  Pig.  2),  Pig.  5  shows  the  period  of  osoillation 
as  a  funotion  of  i  for  synchronous  orbits. 

The  lunar  orbital  plane  is  not  fixed  in  the  ecliptic  but  is  itself 
regressing  around  the  pole  of  the  eoliptio  with  a  period  of  18,6  years  at  an 
inclination  of  5°9'  to  the  eoliptio.  Par  each  position  of  the  lunar  plane  it 
would  be  possible  to  find  the  eigenvalues  X^ ,  Xg,  X^  (X^  ♦  0)  and  the  direotions 

of  the  principal  axes,  but  these  are  varying  all  the  time.  However  it  is  simple 
to  integrate  (16)  numerically,  and  we  have  performed  this  integration  for  an 
initially  equatorial  oiroular  satellite  orbit  for  four  different  initial  posi¬ 
tions  of  the  lunar  plaits.  3inoe  the  polar  coordinates  of  R  with  respeot  to  an 
equatorial  system  are  i  and  (fi-x/2),  where  i  and  (1  are  the" inclination  and  the 
right  asoension  of  the  asoending  node  on  the  equatorial  plane,  these  provide  a 
useful  way  of  illustrating  the  results.  We  have  given  in  Pig.  6  a  polar  plot  of 
the  coordinates  i  and  (n-x/2)  of  the  pole  R.  The  four  different  ourves  refer  to 
the  different  initial  conditions  given  by  ft'  ■  90°  (July  1964),  O'  ■  0°  (Marsh 
1969),  O'  .  270°  (Nov.  1973),  O'  .  180°  (July  1978),  where  O'  is  the  longitude 
of  the  asoending  node  of  the  lunar  ozbit  on  the  plane  of  the  eoliptio  measured 
from  vernal  equinox.  The  symbols  on  the  curves  mark  the  subsequent  positions  of 
the  pole  R  at  intervals  of  5  years  from  the  launch  da to.  In  general  the  tra- 
jeotory  doss  not  quits  pass  through  the  Berth's  axis  after  one  oirouit,  end  the 
time  of  olosost  passage  varies  between  shout  52,5  snd  53.5  years. 
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The  most  significant  feature  revealed  by  Pig.  6  is  that  the  different 
trajectories  interseot  after  about  28  years,  and  again  after  about  56  years. 

This  arises  from  the  near-o  crane  insurability  of  the  period  of  lunar  preoession 
(18,6  years)  and  the  period  of  the  notion  of  the  satellite  cubital  plane 
(~52. 9  years).  It  is  possible  to  give  a  perturbation  treabnent  for  the  dis¬ 
placement  of  the  pole  R  from  the  position  it  would  have  if  the  lunar  orbit  lay 
in  the  eoliptio.  If  the  pole  R  is  not  too  far  away  from  03,  there  are  two  large 
terms  in  the  displacement  with- frequencies  ft1  ±  u,  where  ft'  is  the  rate  of 
regression  of  the  lunar  node  on  the  eoliptio  and  u  is  a  mean  rate  of  preoession 
of  R  around  03.  Numerically 


ft*  i  u  =  (19.34  ±  6.80)°/yr.  >  26.14°/yr.  or  12.54°/yr., 

which  gives  periods  of  13.8  years  and  28,7  years.  Consequently  the  two  large 
terms  in  the  displacement  both  vanish  first  after  about  28  years, 

7  C0NCL0SIQN3 

The  long-term  evolution  of  a  oiroular  satellite  orbit  at  a  distance  of  a 
few  Earth  radii  has  been  investigated.  Only  the  leading  terms  are  retained  in 
the  disturbing  function,  whioh  is  averaged  over  the  mean  anomalies  of  the 
satellite,  the  3un  and  the  Moon.  Negleot  of  the  third  harmonic  in  the  lunar 
disturbing  function  sets  an  upper  limit  to  the  radius  of  the  orbit  for  which 
the  theory  is  valid.  On  the  other  hand  the  period  of  the  predicted  motion  must 
be  appreciably  longer  than  a  year  for  the  double  averaging  procedure  to  be 
justifiable}  this  sets  a  lower  limit  to  the  orbital  radius.  For  those  reasons 
the  analysis  is  most  useful  far  orbits  at  distances  between  3  and  10  Earth 
radii. 


The  motion  ooneists  of  simultaneous  preoession  about  the  Earth's  axis,  the 
pole  of  the  eoliptio  and  the  pole  of  the  lunar  orbit,  and  oan  be  described  by  a 
general  solution  if  those  axes  are  assumed  fixed.  The  motion  of  the  pole  of  the 
orbit  is  analogous  to  the  motion  of  the  instantaneous  axis  of  rotation  of  a 
rigid  body  moving  under  no  forces.  For  a  oiroular  orbit  of  given  else,  there 
are  three  mutually  perpondioular  directions  in  which  the  pole  of  the  orbit  will 
remain  at  rest,  two  of  whioh  are  stable  positions.  In  general  the  projection  of 
the  pole  of  the  orbit  on  the  unit  sphere  follows  a  fixed  tra jeotory  whioh 
enoiroles  one  or  other  of  the  two  stable  positions;  in  faot  the  unit  sphere  is 
divided  into  two  different  regions  bounded  by  a  pair  of  planes  through  the  third 
and  unstable  axis  and  the  tra  jootory  must  lie  entirely  in  one  or  the  other  of 
these  two  regions.  In  practice  the  larger  region  is  the  more  important  and  the 
o corresponding  stable  position  of  the  plane  is  the  well-known  Laplaoian  or  proper 
plane.  The  significance  of  the  proper  plane  is  that  the  pole  of  an  orbit 
inclined  at  a  small  angle  to  it  will  regress  around  the  pole  of  the  proper  plane 
at  nearly  oomtant  rate  and  inclination. 

The  analysis  has  been  applied  to  Earth  satellites  assisting  the  lunar  orbit 
lies  in  the  eoliptio.  The  inclination  of  the .proper  plane  to  the  equator 
increases  from  0.19°  at  3  Earth  radii  to  I8.8°at  10  Earth  radii,  and  the  semi¬ 
angle  between  the  bounding  planes  reaches  a  maximum  value  of  nearly  12°  at  about 
7.7  Earth  radii,  as  shown  in  Fig.  2.  As  the  tra  jeotory  of  the  pole  of  the  oxbit 
shrinks  towards  either  of  the  stable  positions  the  period  deoreases  to  one  or 
other  of  two  limiting  values.  For  a  given  size  of  orbit,  the  limiting  value  at 
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the  pole  of  the  proper  plane  is  the  shortest  possible  period)  this  limiting 
period  attains  a  marimm  value  of  about  70  years  near  8. 9  Earth  radii  as  shown 
in  Pig.  4*  The  limiting  period  around  the  other  stable  stationary  position 
inoreases  nearly  linearly  with  radial  distance  from  121  to  403  years  between 
3  and  10  Earth  radii. 

Synchronous  orbits  havo  been  oonsiderod  in  more  detail.  The  period  of  the 
long-term  motion  is  given  for  all  initial  positions  in  Pig. 5,  and  the  limiting 
values  of  the  period  are  32.5  and  267  years.  The  effoot  of  the  motion  of  the 
lunar  orbit  on  an  initially  equatorial  orbit  at  synchronous  height  has  been 
found  by  numerical  integration  for  four  differont  initial  positions  of  the  lunar 
orbit,  i,e.  four  difforont  launch  dates,  and  is  illustratod  in  Pig. 6. 
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